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We present a study of the average transverse momentum (pt) fluctuations and p t correlations 
for charged particles produced in Cu+Cu collisions at midrapidity for ^/snn = 62.4 and 200 GeV. 
These results are compared with those published for Au+Au collisions at same energies, to explore 
the system size dependence. In addition to the collision energy and system size dependence, the 
Pt correlations results have been studied as functions of the collision centralities, the ranges in p±, 
the pseudo-rapidity r\, and the azimuthal angle (j>, for which the correlations are measured. The 
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square root of the measured pt correlations when scaled by mean-p t are found to be independent 
of both colliding beam energy and system size studied. The transport based model calculations are 
found to have a better quantitative agreement with the measurements compared to models which 
incorporate only jet-like correlations. 

PACS numbers: 25.75. Gz 



I. INTRODUCTION 

The study of event-by-event fluctuations and correla- 
tions is an important tool to understand thermalization 
and phase transitions in heavy-ion collisions [3-Q. Non- 
monotonic change in transverse momentum (p t ) corre- 
lations as a function of ccntrality and/or the incident 
energy has been proposed to be a possible signal of the 
Quark Gluon Plasma (QGP) formation The QGP is 
believed to be formed at the early stage of high energy 
heavy-ion collisions when the system is hot and dense. 
As the time passes, the system dilutes, cools down and 
hadronizes. Fluctuations are supposed to be sensitive 
to the dynamics of the system, especially at the QGP 
to hadron phase transition 0, H|. Alternatively, analy- 
ses at the Relativistic Heavy Ion Collider (RHIC) based 
on p t auto-correlations (the inversion of the mean trans- 
verse momentum fluctuations in different pseudorapidity 
and azimuthal angle difference regions of the produced 
particles) indicate that the basic correlation mechanism 
could be dominated by the process of parton fragmen- 
tation 0]. Thus fluctuation measurements are proposed 
to be an important tool in understanding the nucleus- 
nucleus collisions 0, H, H, I8l4l2| . In addition, under the 
assumption that thermodynamics is applicable to heavy- 
ion collisions, fluctuations in various observables could 
be related to thermodynamic properties of the matter 
formed. For example, the event-by-event (pt) could be 
related to the temperature fluctuations [lj-[l7j]. 

The study of event-by-event fluctuations of various 
quantities such as event- wise mean transverse momentum 
((pt)), charged track multiplicity, and conserved quanti- 
ties such as net-baryon and/or net-charge are considered 
to be some of the main probes for the search of the critical 
point in the QCD phase diagram [r7l - l26| . One expects 
enhanced fluctuations in the above observables when the 
system passes through the vicinity of the critical point. 
Recent results from the SPS experiments show that the 
energy dependence of transverse momentum fluctuations 
do not show the increase expected for freeze-out near the 
critical point of QCD (2tJ. However, when these fluctu- 
ations arc studied as a function of the system size (col- 
liding C+C, Si+Si, Pb+Pb), enhanced fluctuations are 
observed in smaller colliding systems [28[. The results 
from the RHIC Beam Energy Scan (BES) [H data for 
the above observables will provide further insights. 

The results presented here are from Cu+Cu collisions 
at y/sNN = 62.4 and 200 GeV using the Solenoidal 
Tracker At RHIC (STAR) and are compared with the 
published results from Au+Au collisions at the same en- 
ergies [30]. This paper describes a systematic study of 



the system size dependence of the transverse momentum 
correlations observed at RHIC. 

The paper is organized as follows. In section II, we 
discuss the detectors, the data set, and the centrality se- 
lection used in the analysis. In section III, we discuss (p t ) 
fluctuations extracted from the (pt) distributions, which 
are compared with mixed events and Gamma distribu- 
tions. Dynamical fluctuations are extracted and pre- 
sented for Au+Au and Cu+Cu collisions at y / s/vw = 62.4 
and 200 GeV. In section IV we discuss the pt correla- 
tions, the dynamical correlations when scaled by (-/Vp art ) 
and (p t ) to understand the centrality and energy depen- 
dence. We also compare the experimental data to various 
model calculations. Finally, we present a detailed study 
of i], (p and pt dependence of the correlations. In sec- 
tion V, we discuss the systematic errors associated with 
the analysis presented. We conclude with a summary in 
section VI. 



II. EXPERIMENT AND DATA ANALYSIS 

The Cu+Cu data were taken using the STAR detec- 
tor with a minimum bias trigger. For the data taken at 
y/sNN = 200 GeV this was done by requiring a coinci- 
dence from the two Zero Degree Calorimeters (ZDCs). 
For the y / S7v r /v = 62.4 GeV data the ZDC is less efficient, 
so a Beam Beam Counter (BBC) coincidence was also 
required. More details about the trigger detectors can 
be obtained from Ref. [3lj . 

The main detector used in this analysis is the Time 
Projection Chamber (TPC) [HJ, which is the primary 
tracking device in STAR. The TPC is 4.2 m long and 4 m 
in diameter and its acceptance spans about ±1.0 units of 
pseudorapidity and full azimuthal coverage. The sensi- 
tive volume of the TPC contains P10 gas (10% methane, 
90% argon) regulated at 2 mbar above atmospheric pres- 
sure. The TPC data are used to determine particle tra- 
jectories, momenta, and particle-type through ionization 
energy loss (dE/dx) (33j . 

The primary vertex of events used in this analysis was 
required to be within ± 30 cm of the geometric center 
of the TPC along the beam axis. This selection process 
yielded about 7.5 M and 15 M minimum bias events for 
Cu+Cu collisions at ^snn = 62.4 and 200 GeV, respec- 
tively. 

The collision centralities are represented by the frac- 
tions of events in the minimum bias inelastic cross sec- 
tion in a collision. In data, the collision centrality is 
determined by using the uncorrected charged track mul- 
tiplicity (A^ PC ), measured in the TPC within \r]\ < 0.5. 
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The various centrality bins are calculated as a fraction of 
this multiplicity distribution starting at the highest mul- 
tiplicities. The centrality classes for Cu+Cu collisions 
at y/sj^ = 62.4 and 200 GeV are 0-10% (most central), 
10-20%, 20-30%, 30-40% 40-50% and 50-60% (most pe- 
ripheral). Each centrality bin is associated with an aver- 
age number of participating nucleons ((iVpart)) obtained 
using Glauber Monte Carlo simulations [34J employing 
the Woods-Saxon distribution for the nucleons inside the 
Cu nucleus. The systematic uncertainties are determined 
by varying the Woods-Saxon parameters and the uncer- 
tainties associated with nucleon-nucleon cross sections. 
Table U lists the N]^ c and (N paTt } values for each cen- 
trality in Cu+Cu collisions at ^/snn = 62.4 and 200 
GeV. Corresponding values for Au+Au collisions can be 
found in Ref. 

In order to have uniform detector performance, a pscu- 
dorapidity cut of < 1.0 is applied to tracks. To avoid 
the admixture of tracks from a secondary vertex, a re- 
quirement is placed on the distance of closest approach 
(DCA) between each track and the event vertex. The 
charged particle tracks are required to have originated 
within one cm of the measured event vertex. The mul- 
tiple counting of split tracks (more than one track re- 
constructed from the original single track) is avoided by 
applying a condition on the number of track fit points 
(A^Fit) used in the reconstruction of the track. We require 
that each included track had a minimum number of 21 
TPC points along the track. The transverse momentum 
range selected for the analysis is 0.15-2.0 GeV/c. 

The errors shown in the figures have the statistical and 
systematic errors added in quadrature unless otherwise 
stated. The statistical errors are small and are within 
symbol sizes. The final systematic errors are obtained 
as quadrature sums of systematic errors from different 
sources as discussed in section V. 



III. (p t ) FLUCTUATIONS 

The p t fluctuations in high energy collisions can be 
measured using the distribution of the event-wise mean 
transverse momentum defined as 



(Pt 



1 N 



(i) 



i=i 



where N is the multiplicity of accepted tracks from the 
primary vertex in a given event and p tt i is the transverse 
momentum of the i th track. The mean-p t distribution is 
compared to the corresponding distribution obtained for 
"mixed events" j3(j. Mixed events for a given central- 
ity are constructed by randomly selecting one track from 
each event in that centrality while keeping the same event 
vertex. The total multiplicity for a given mixed event is 
kept the same as in the real data. 

Figure [TJ shows the event- by- event mean-p t distribu- 
tions for 0-10% Cu+Cu collisions at (a) -y/sjvw = 200 and 



TABLE I: The Nj? c values and average number of partici- 
pating nucleons ((A pa rt)) for various collision centralities in 
Cu+Cu collisions at ^/snn = 62.4 and 200 GeV. 



% cross 
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FIG. 1: (color online) Event-by-event (p t ) distributions for 
data and mixed events in central Cu+Cu collisions at (a) 
y/SNN — 200 and (b) 62.4 GeV. The curves (solid for data 
and dotted for mixed events) represent the V distributions. 
The errors shown are statistical. 



(b) 62.4 GeV. The solid symbols represent the (pt) distri- 
butions for data; the open symbols represent (pt) distri- 
butions for mixed events. The distributions are similar 
for other centralities. The mixed events provide reference 
to statistical fluctuations in the data. Any fluctuations 
observed in data beyond these statistical fluctuations are 
referred to as non-statistical or dynamical fluctuations in 
this paper. For both data and mixed events, as we go 
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TABLE II: Gamma distribution parameters for event-by- 
event mean-(p t } distributions for data and mixed events in 
central (0-10%) Cu+Cu collisions at ^/sWn = 62.4 and 200 
GeV. 



Collision 


Q 


P (xl0 _a 




a 


type (AA) 




GeV/c) 


(GeV/c) 


(GeV/c) 


Cu 200 (Data) 


476 


1.22 


0.5805 


0.02660 


Cu 200 (Mixed) 


634 


0.92 


0.5807 


0.02310 


Cu 62.4 (Data) 


358 


1.56 


0.5603 


0.02960 


Cu 62.4 (Mixed) 


457 


1.23 


0.5602 


0.02621 



from central to peripheral collisions, the mean (fx) of the 
distributions decreases while the width (a), measured as 
the root mean square, increases. Moreover, it is seen that 
the (jpt) distributions for data are wider than those for 
mixed events, suggesting the presence of non-statistical 
fluctuations in Cu+Cu data for both 62.4 and 200 GeV 
collisions. 

The curves in Fig. [T] represent the Gamma (T) distri- 
butions for data (solid lines) and mixed events (dotted 
lines). The T distribution [3, [35| is given by 



/(*) 



jM—lg—tc/P 

T(a)(3 a ' 



(2) 



where x represents the (p t ); a = /i 2 /c 2 and /3 = cr 2 //x. 

According to Ref. [HI the V distribution is one of the 
standard representations of the inclusive single particle 
Pt distribution. The author also suggests that the quan- 
tity a/ (N c b), should be ~ 2, and the quantity (3 x {Nek) 
representing the inverse slope parameter may be referred 
to as the temperature of the system. Here (N c h) is the 
average charged particle multiplicity in a given centrality 
bin. We find that a/(N ch ) for Cu+Cu 0-10% central col- 
lisions at 200 GeV is 2.04 and that at 62.4 GeV is 2.18. 
The respective (3 x (N c h) values are 0.284 GeV/c and 
0.256 GeV/c. The a/{N ch ) and /3 x (N ch ) for 0-5% cen- 
tral Au+Au collisions at 200 GeV were found to be 1.93 
and 0.299 GeV/c H3. TablelQlists Gamma distribution 
parameters for event-by-event (pt) distributions for data 
and mixed events in central (0-10%) Cu+Cu collisions 
at yi^" = 62.4 and 200 GeV. The values for Au+Au 
collisions can be found in Ref. [30j |. For Cu+Cu colli- 
sions at 200 GeV, a/{N c k) varies from 2.04 to 2.11 from 
central to peripheral collisions. However, for Cu+Cu at 
62.4 GeV, its values varies from 2.18 to 2.27 from central 
to peripheral collisions. The errors on values of a and /3 
are of the order of 13-18% and 9-12%, respectively, for 
Cu+Cu collisions. 

The non-statistical or dynamical fluctuations in mean- 
Pt are quantified using a variable <7d y n [IB, defined 
as 



0" dyn 



"date 



(3) 



where //data and /x m i x are the means of the event-by-event 



(pt) distributions for data and mixed events, respectively. 
Similarly, rjdata and cr mix are respectively the root mean 
square deviations in (p t ) distributions for data and mixed 
events. Figure [5] shows the dynamical fluctuations (cdyn) 
in mean-p t plotted as a function of (V par t). The results 
are shown for Cu+Cu collisions at ^/sjvw = 62.4 and 200 
GeV, and are compared with the results from Au+Au 
collisions at ^/s^n = 62.4 and 200 GeV. The dynami- 
cal (pt) fluctuations are similar in Au+Au and Cu+Cu 
collisions at similar values of (Af part ). The fluctuations 
decrease as the (N part ) increases. The dynamical fluctu- 
ations are also independent of the collision energy and 
are found to vary from ~ 2% to ~ 5% for (N palt ) less 
than ~+20. For (N pai t) greater than ~150, the dynami- 
cal fluctuations are smaller and vary from 1% to 2.5%. 
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FIG. 2: (color online) Comparison of dynamical (p t ) fluctu- 
ations in Au+Au and Cu+Cu collisions at ^/snn = 62.4 and 
200 GeV as a function of number of participating nucleons. 



IV. Pt CORRELATIONS 

The non-statistical or dynamical fluctuations can also 
be analyzed by using the two-particle transverse momen- 
tum correlations [30j. The two-particle pt correlations 
are studied using the following equation [30(: 



(Ap M Ap 



1 



t,j! 



E 



Advent f-i N k (N k -l)' 



(4) 



where C k is the two-particle transverse momentum co- 
variance for the k th event, 

ft = E E (PM-<<Pt»)(PW-<<Pt»), (5) 

i=l j=X,ijtj 

Pt y t is the transverse momentum of the i track in the 
k event, Nk is the number of tracks in the k th event, 
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and iV event is the number of events. The overall event 
average transverse momentum (((pt))) is given by 

/Ncvcnt \ 

«P*» = ( E <P*>* J /^-ent, (6) 

where (pt) fe is the average transverse momentum in the 
k th event. It may be noted that for mixed events, there 
will be no dynamical fluctuations and the pt correlations 
will be zero. Equation (|4]) is used to obtain the pt correla- 
tions in Cu+Cu collisions at ^/snn = 62.4 and 200 GeV. 
These results are compared with the published results 
from Au+Au collisions at similar energies (30| to investi- 
gate the system-size and collision energy dependence of 
the pt correlations in heavy- ion collisions at RHIC. 

The p t correlation values may be affected by the de- 
pendence of the correlations on the size of the centrality 
bin and by changes induced due to a change in ((p*)) . We 
removed this dependence by calculating ((pt)) as a func- 
tion of (A c h). We fitted this dependence using a suitable 
polynomial function and used this fit in Eq. ([5]). All re- 
sults presented in this paper have been corrected for this 
effect. 

Figure [3J (a) shows the pt correlations plotted as func- 
tion of (iVpart) for Cu+Cu and Au+Au collisions at 
^/Jnn = 62.4 and 200 GeV. We observe pt correlations 
that decrease with increasing (N part ) for Au+Au and 
Cu+Cu at both energies. The decrease in correlations 
with increasing participating nucleons could suggest that 
correlations are dominated by pairs of particles that orig- 
inate from the same nucleon-nucleon collisions, and they 
get diluted when the number of participating nucleons 
increases [3(3]. 

A. Scaled pt Correlations 

The decrease in pt correlations with increasing (V par t) 
seen in Fig. [3J (a) may be related to a system volume 
dependence characterized by (iVpart). To remove this 
volume dependence we multiply the pt correlations by 
(Af p art)/2 as shown in Fig. [3J (b). The results are shown 
for Cu+Cu and Au+Au collisions at ^/snn = 62.4 and 
200 GeV. It is observed that this measure of pt correla- 
tions increases quickly with increasing collision centrality 
for both Cu+Cu and Au+Au collisions and saturates for 
central Au+Au collisions. The saturation of this quan- 
tity might indicate effects such as the onset of thermal- 
ization [2^|, the onset of jet quenching 0, [Tl|, or the 
saturation of transverse flow in central collisions fo8j . It 
seems that for Cu+Cu collisions this quantity is larger 
than for Au+Au collisions which might indicate more 
correlations for the smaller systems. But since the errors 
are large it is hard to make this conclusion. 

The correlation measure {Ap t ^ Ap t j) ma y change due 
to changes in ((pt)) with incident energy and/or colli- 
sion centrality. To address these changes, we studied 
the square root of the measured correlations scaled by 




5 r (d) i 

Ofri = 

50 100 150 200 250 300 350 

<N Part > 



FIG. 3: (color online) (a) p t correlations, (b) p± correlations 
multiplied by (A r par t}/2, (c) square root of p t correlations 
scaled by {{pt}), and (d) square root of pt correlations multi- 
plied by (iVpart)/2 and scaled by {{pt)), plotted vs. {N palt )- 
Results are compared between Cu+Cu and Au+Au collisions 
at ^/snn = 62.4 and 200 GeV. Au+Au data has been taken 
from the Ref. [U . 



{{pt}}- Figure [3j (c) shows the corresponding quantity 
yj {Apt.i Ap t j) I ((Pt)) plotted as a function of collision 
centrality for Cu+Cu and Au+Au collisions at ^/snn = 
62.4 and 200 GeV. It is observed that the correlation 
scaled by ((pt)) is independent of collision system size and 
energy, but decreases with increasing collision central- 
ity. The combined effect of multiplying p t correlations by 
(iVpart) /2 and scaling with ((pt)) is shown in Fig.|3j(d). It 
seems that this quantity ^(iY part /2) (Ap M Ap t j)/{{p t )) 
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increases with (-ZV part ) and shows saturation for central 
Au+Au collisions, but is independent of collision system 
and energy. 



B. Model Comparisons 

It is interesting to compare the above results with the- 
oretical model calculations to understand the physical 
mechanism behind these measurements. We compared 
our results with some widely used models in heavy-ion 
collisions such as Ultra Relativistic Quantum Molecular 
Dynamics model - UrQMD [H], A Multi Phase Trans- 
port model - AMPT (default and string-melting) jio| . 
and Heavy- Ion Jet Interaction Generator - HIJING (with 
jet quenching switched off and on) [4l|. We have used 
UrQMD ver.2.3, AMPT ver.l.II for default, and ver.2.11 
for AMPT string-melting. 

HIJING is a pcrturbative QCD-inspired model that 
produces multiple minijct partons; these later get trans- 
formed into string configurations and then fragment 
to hadrons. The fragmentation is based on the Lund 
jet fragmentation model [Hj]. A parametrized parton- 
distribution function inside a nucleus is used to take into 
account parton shadowing. 

The AMPT uses the same initial conditions as in HI- 
JING. However, the minijet partons are made to un- 
dergo scattering before they are allowed to fragment 
into hadrons. The string-melting (SM) version of the 
AMPT model (labeled here as AMPT Melting) is based 
on the idea that for energy densities beyond a critical 
value of 1 GeV/fm 3 , the system should be de-confined 
and strings (or hadrons) decomposed into their partonic 
components. This is done by converting the mesons to 
a quark-antiquark pair, baryons to three quarks, and so 
on. The scattering of the quarks is based on a parton 
cascade. Once the interactions stop, the partons then 
hadronize through the mechanism of parton coalescence. 
The UrQMD model is based on a microscopic transport 
theory where the phase-space description of the reactions 
is important. It allows for the propagation of all hadrons 
on classical trajectories in combination with stochastic 
binary scattering, color string formation, and resonance 
decay. 

Figure 2] shows the comparison of p t correlations be- 
tween data ((a) and (c) for Au+Au 200 GeV, and (b) 
and (d) for Cu+Cu 200 GeV) and the models described 
above. We observe that the transport based UrQMD 
model calculations are found to have a better quanti- 
tative agreement with the measurements compared to 
models which incorporate only jet-like correlations as in 
HIJING. HIJING gives similar dependence on (Apart) 
but under-predicts the magnitude. Inclusion of the jet- 
quenching effect in the HIJING leads to a smaller value 
of the correlations in central collisions. A multi-phase 
transport model calculation incorporating coalescence as 
a mechanism of particle production also compares well 
with data for central collisions. However the default 
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FIG. 4: (color online) Comparison of pt correlations between 
data and models. Au+Au data has been taken from the 
Ref. [3(j. The curves represent different model calculations. 
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FIG. 5: (color online) p t correlations for varying rapid- 
ity acceptance (\r)\ < 0.25,0.5,1.0) for Cu+Cu collisions at 
V^jviv = 200 and 62.4 GeV. 
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FIG. 6: (color online) pt correlations for varying azimuthal 
acceptance for Cu+Cu collisions at ^/snn = 200 (panels (a) 
and (c)) and 62.4 GeV (panels (b) and (d)). 
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version of this model which incorporates additional ini- FIG. 7: (color online) p t correlations for different p t ranges 
tial and final state scattering relative to HIJING yields shown for (a) Cu+Cu 200 GeV and (b) Cu+Cu 62.4 GeV. 
smaller values of correlations. 



C. rj and (f> Dependence 

The rj and </>- dependences of p t correlations were also 
studied. Figure [5] shows the pt correlations plotted as 
a function of increasing rapidity acceptance for Cu+Cu 
collisions at (a) ^/sWn = 200 and (b) 62.4 GeV. The 
data points for ccntralitics 0-10%, 20-30%, and 40-50% 
arc shifted by 0.01 in A?//2 for clarity. The correlations 
are almost independent of the A77 window for most cen- 
tral collisions. For peripheral collisions, the correlations 
show a slight rapidity dependence with maximum value 
for -0.25 < 77 < 0.25. 

Figure [S] shows the pt correlations for varying az- 
imuthal angle windows for Cu+Cu collisions at (a) 
y/sW^ = 200 and (b) 62.4 GeV. The data points for 
4> windows A</> : 30° and 90° are shifted by 0.5 units 
in -/Vp ar t for clarity. We observe that the correlations are 
highest for the Ac/) = 15° (among the cases studied). The 
correlations decrease with increasing collision centrality 
for a given Acj) window as expected. The correlations 
multiplied by (jV par t)/2 (Fig. 1(c) for Cu+Cu 200 GeV 
and Fig. [5] (d) for Cu+Cu 62.4 GeV) seem to increase 
and then saturate with increasing collision centrality. 



D. pt Dependence 

Figure [7] shows the correlations as a function of col- 
lision centrality for different p t regions for Cu+Cu col- 
lisions at (a) ^/IWn = 200 and (b) 62.4 GeV. The dif- 
ferent pt ranges used are shown. The selection of these 
Pt ranges is done in order to see the dependence of p t 
correlations among the lower p t , higher p t , and across 
higher and lower p t tracks. The data points for p t ranges 
(in GeV/c): 0.15-0.5, 0.5-2.0, 0.15-0.35, and 0.35-2.0, are 
shifted by one unit in (-/V part )- 

We observe that pt correlation is maximum (minimum) 
for charged particles whose p t is between 0.15-2.0 GeV/c 
(0.15-0.35 GeV/c). The p t correlation values are small 
and fairly independent of p t if a lower p t bound for the 
particles studied is fixed at 0.15 GeV/c and the upper 
Pt bound is progressively increased up to 0.50 GeV/c. 
When the analysis is carried out by keeping the higher 
Pt bound fixed at 2.0 GeV/c and subsequently decreasing 
the lower p t bound to 0.15 GeV/c, the correlation values 
are found to increase. 

Figure [H] shows the variation of p t correlations as a 
function of (N part ) for different p t windows as calcu- 
lated using the (a) AMPT (string-melting), (b) UrQMD, 
and (c) HIJING (no jet quenching) model calculations 
for Cu+Cu collisions at 200 GeV. The AMPT calcula- 
tions show pt correlations that are similar to those ob- 
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FIG. 8: (color online) Correlations for varying p t ranges for 
different model calculations in Cu+Cu 200 GeV: (a) AMPT 
(SM), (b) UrQMD, and (c) Hijing (No JQ). 



served in data for corresponding variations in the pt win- 
dows. The trend of the correlation values shown by both 
UrQMD and HIJING is similar to what is seen in the 
data for the low-p t windows where the lower bound is 
fixed at 0.15 GeV/c and the upper bound is increased 
from 0.35 GeV/c to 0.50 GeV/c. However, for UrQMD, 
if the higher p t bound is fixed at 2.0 GeV/c and the 
lower p t bound is subsequently decreased to 0.15 GeV/c, 
Pt correlation values remain similar. For the same case, 
HIJING shows decrease in p t correlation values when the 
lower p t bound is decreased to 0.15 GeV/c. This is just 
opposite of what was observed in data as seen in Fig. [7J 

Because we calculate correlations for different p t 
ranges, we have examined the pt acceptance effects on 
the observed pt correlations. The correlation values in 
different pt ranges may depend on the pt range size and 
on the fluctuations in the pt spectrum slope in that pt 
range. It is, therefore, important to see the effect of 
slope fluctuations on the correlation values in different p t 
ranges. The pt correlations can be formulated in terms of 
the fluctuations in the inverse slope parameter (effective 
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FIG. 9: (color online) Top panel: Uncorrected inclusive 
charged particle pt spectrum for 0-10% collision centrality 
for Cu+Cu 200 GeV (open circles). The distribution is fit 
to the exponential function: Ae~ Pt ^ T (red dashed-line). Er- 
rors are statistical. Bottom panel: Correlations scaled by 
(d(p t )/dT) 2 T 2 vs. (A^part) for different p t ranges in Cu+Cu 
200 GeV. 



temperature) by the following relation [11 



(Apt Apt 



d(j>t) 



-i 2 



dT 



AT 2 , 



(7) 



where AT 2 describes the fluctuation in the inverse slope 
parameter. The dependence (p t (T)) can be obtained 
from the function that describes the inclusive uncorrected 
pt spectrum for the desired pt range in the following man- 
ner. Figure |H] (top panel) shows the measured inclusive 
Pt spectrum for 0-10% collision centrality in Cu+Cu col- 
lisions at 200 GeV. The dotted curve represents the ex- 
ponential fit of the form F(p t ) = Ae~ Pt / T , that is fit 
to these measurements for the pt range 0.35-2.0 GeV/c. 
The expression for (p t ) and (d(pt)/dT) can be obtained 
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using the following relation: 



Pt F (Pt)dpt 



p t F(p t )dpt 



(8) 



which gives, 



(Pt) =2T + 



a 2 e -a/T _ b 2 e -b/T 

(a + T)e- a / T - {b + T)e- b ' T ' 



(9) 



Here, a and b are respectively the lower and upper limits 
of a given p t range. The derivative of (pt) with respect 
to T is obtained as: 



dipt) 
dT 

where 
A = 



Ae -(a+b)/T + a 2 e -2a/T + 

2 ~ [( a + r)e- a / T -(6 + T)e- b /r]2 ' ( 10 ) 



a6(6-a) 2 (6 2 -a 2 )(6-a) 



y2 



(a 2 +6 2 ). (11) 



Using the lower and upper limits of a given p t range, the 
corresponding T from the spectrum fit for each pt range 
and collision centrality, we can obtain (d(p t )/dT) 2 for 
every p t range in each centrality. 

Figure [9] (bottom panel) shows the measured two- 
particle p t correlations scaled by (d(p t ) / dT) 2 T 2 vs. 
(iVp ar t) for different p t ranges in Cu+Cu 200 GeV. The 
data points for p t ranges (in GcV/c): 0.4-0.8, 0.5-2.0, 1.2- 
1.6, and 0.35-2.0, arc shifted by one unit in (-/Vp art ) for 
clarity. We observe that scaled pt correlations for differ- 
ent pt ranges become similar and show little dependence 
on the collision centrality. This study seems to suggest 
that the difference in the p t correlations for different p t 
ranges is due to the fluctuations in slope of the p t spec- 
trum in those pt ranges. 



V. SYSTEMATIC UNCERTAINTIES 

Systematic errors on the mean (fj,) and root mean 
square deviations (a) in the (j>t) distributions (discussed 
in section III), and p t correlations (discussed in sec- 
tion IV) are mainly evaluated by varying the different 
cuts used in the analysis, re-doing analysis using these 
changed cuts and determining the resulting changes in 
the values of fi, a, and the p t correlations. The difference 
is taken as the systematic error due to a particular analy- 
sis cut. The resulting systematic uncertainties, described 
below, are shown in Tables IIIII and IIVI as a percentage 
of the result (p, a, and the pt correlations) for various 
centralities for Cu+Cu collisions at both 62.4 and 200 
GeV. 

To study the effect of the z- vertex (V z ) cut we increased 
the V z acceptance to ± 50 cm from the default value of 
±30 cm. We observed no change in \x or a or in the pt 
correlations when using the wider V z . 



TABLE III: Systematic errors on /i and a in event-wise {pt) 
distributions as described in section V for various collision 
centralities in Cu+Cu collisions at ^/snn = 62.4 and 200 
GeV. 





Cu+Cu 62.4 GeV 


Cu+C 


u 200 GeV 




% cross 


DCA 


Afrit 


DCA 


Afrit 


DCA 


Ar Fit 


DCA 


Afrit 


section 


M(%) 


M(%) 


a{%) 


a(%) 


M(%) 


M(%) 


a(%) 


a(%) 


0-10 


3.6 


0.4 


7.9 


1.0 


3.4 


0.24 


5.7 


1.2 


10-20 


3.6 


0.4 


7.1 


1.0 


3.3 


0.23 


5.2 


1.2 


20-30 


3.6 


0.4 


6.6 


1.0 


3.3 


0.23 


5.2 


1.1 


30-40 


3.6 


0.4 


6.2 


1.0 


3.3 


0.20 


4.7 


1.2 


40-50 


3.6 


0.4 


6.0 


1.0 


3.2 


0.20 


4.9 


1.1 


50-60 


3.6 


0.4 


6.0 


1.0 


3.2 


0.22 


4.4 


1.2 



TABLE IV: Systematic errors on pt correlations as described 
in section V for various collision centralities in Cu+Cu colli- 
sions at yfsWN = 62.4 GeV and 200 GeV. 





Cu+Cu 62.4 GeV 


Cu+Cu 200 GeV 


% cross 
section 


DCA 

(%) 


Afr it 
(%) 


Poly. 

(%) 


Low pt 

(%) 


DCA 

(%) 


Afr t 
(%) 


Poly. 

(%) 


Low pt 

(%) 


0-10 


30 


1.4 


1.9 


7.2 


16 


0.05 


1.9 


22 


10-20 


23 


0.8 


3.6 


13.1 


13 


0.09 


3.6 


3.2 


20-30 


18 


0.6 


1.9 


3.4 


13 


1.1 


1.9 


12.3 


30-40 


17 


1.0 


0.004 


9.0 


8 


1.2 


0.004 


9.7 


40-50 


19 


3.0 


0.009 


1.0 


10 


2.0 


0.009 


8.4 


50-60 


20 


3.6 


0.009 


4.0 


7 


5.0 


0.009 


8.3 



The effect of the cuts used to suppress background 
tracks was studied by changing the DCA cut from the 
default, DCA < 1 cm, to DCA < 1.5 cm and sepa- 
rately, changing the required number of fit points along 
the track, AT Fit , from the default V Fit > 20 to V Fit > 
15. The resulting systematic errors due to these changes 
are listed in Tables IIIII and IIVI in the columns labeled 
"DCA" and "AW ■ 

The effect of the size of the centrality bin on the p t cor- 
relations was addressed by fitting ((pt)) as a function of 
(V c h) (see section IV). To determine the systematic un- 
certainty associated with this process we used different 
polynomial functions to fit ((pt)) vs. (V c h). The sys- 
tematic errors associated with this correction are listed 
in Table UVl in the columns labeled "Poly". 

We estimated the systematic uncertainty on the pt cor- 
relations that may be associated with the application of 
the low-pi cut by removing this pt cut in the HIJING[4l| 
model calculations. The estimated systematic errors are 
shown in Table [TVl in the columns labeled "Lowpt". 

The p t correlations may also include the short range 
correlations such as Coulomb interactions and Hanbury 
Brown- Twiss (HBT) correlations. These correlations 
usually dominate between the pairs of particles having 
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relative transverse momentum less than 100 MeV/c. The 
effect of these short-range correlations on the measured 
p t correlations is seen by calculating p t correlations after 
removing the pairs of particles with relative momentum 
(Pi — Pj)> l ess than 100 McV/c. We observed that the p t 
correlations are reduced by maximum of 6% when short- 
range correlations are excluded. 

The p t correlations also include the contributions from 
resonance decays and charge ordering. These correla- 
tions are obtained for pairs of particles having like (H — h 
or ) and unlike (-) — ) charges with respect to inclu- 
sive charged particles. A maximum of 15% decrease in 
the correlations is observed for pairs of particles with 
like charges and about 12% increase is observed for pairs 
with unlike charges with respect to the correlations for 
inclusive charged particle pairs for Cu+Cu collisions at 
^/sWn = 62.4 GcV and 200 GeV. 

VI. SUMMARY 

We have reported the measurements of pt fluctuations 
in Cu+Cu collisions in the STAR detector at RHIC for 
y/sjVTV = 62.4 and 200 GeV, and compared with Au+Au 
collisions at the same energies to investigate the sys- 
tem size dependence. We have observed dynamical pt 
fluctuations by comparing data to mixed events in both 
Cu+Cu and Au+Au collisions at these two beam ener- 
gies. Moreover, for similar mean number of participating 
nucleons, the pt fluctuations are observed to be compara- 
ble for the Cu+Cu and Au+Au systems, suggesting that 
the system size has little effect. In addition, we have 
compared pt correlation measurements for Cu+Cu colli- 
sions at ^snn = 62.4 and 200 GeV with the published 
Au+Au measurements. We found that for both Cu+Cu 
and Au+Au systems the p t correlation decreases with in- 
creasing (iVp art ) at both beam energies. The dilution of 
the p t correlations with increasing (A par t) could be due to 
the decrease in correlations that are dominated by pairs 
of particles that originate from the same nucleon-nucleon 
collision. A similar observation is found with other fluc- 
tuation and correlation observables like K/tt [37J and net- 
charge [43j ] . 

The pt correlations were found to scale with number 
of participating nucleon pairs for (A par t) > 100 (when 
the system size is larger than that of central Cu+Cu 
collisions). This might indicate the onset of thermal- 
ization [26|, the onset of jet quenching 0, [ll|, or the 
saturation of transverse flow in central collisions [38| . 
The square root of the pt correlations normalized by 
event-wise average transverse momentum for Cu+Cu 
and Au+Au collisions is similar for systems with similar 
(Ap ar t) and is independent of the beam energies studied. 

The results described above are compared to predic- 
tions from several relevant model calculations. The 



transport based UrQMD model calculations are found to 
have a better quantitative agreement with the measure- 
ments compared to models which incorporate only jet- 
like correlations as in HIJING. HIJING gives similar de- 
pendence on (Apart), but under-predicts the magnitude. 
Inclusion of the jet-quenching effect in the HIJING leads 
to a smaller value of the correlations in central collisions. 
A multi-phase transport model calculation incorporating 
coalescence as a mechanism of particle production also 
compares well with data for central collisions. 

When studying the p t correlations for different p t in- 
tervals, the correlations appear to be small and fairly 
independent of pt interval, if we fix the lower p t bound 
at 0.15 GeV/c and progressively increase the higher p t 
bound up to 0.50 GeV/c. This suggests that correlations 
are weak for low-pt particles. This \ow-pt trend observed 
in the data is also seen in UrQMD, AMPT, and HIJING 
models. When the analysis is carried out keeping the 
higher p t bound fixed at 2.0 GeV/c and progressively 
decreasing the lower p t bound to p t = 0.15 GeV/c, the 
correlation values in data are found to increase. This 
suggests that high p t particles are more correlated with 
low-p t particles. The AMPT model shows a rather sim- 
ilar variation of pt correlations for different p t intervals, 
as observed in data. The UrQMD model calculations, 
however, show no such variations in correlations for the 
different pt intervals with higher p t bound fixed at 2.0 
GeV/c. Finally, we note that the HIJING model calcula- 
tions give pt correlations that decrease with a decrease in 
the lower p t bound for intervals with fixed higher p t (= 
2.0 GeV/c) bound that is opposite to what is observed 
in data. We tried to understand the changes in p t corre- 
lations in different p t intervals as observed in data. We 
observed that the resultant fluctuations after considering 
event- by- event variation in the slope of the pt spectra for 
different p t bins were all of similar order. 

We have also shown the variation of the pt correlation 
with the change in the accepted range of pscudorapidity 
and azimuthal angle of the produced particles. We ob- 
served that the correlation values increase when the 77- 
and ^-acceptance are reduced. 
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